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We investigate the nonlinear interaction between a relativistically strong laser beam and a plasma 
in the quantum regime. The collective behavior of the electrons is modeled by a Klein-Gordon equa- 
tion, which is nonlinearly coupled with the electromagnetic wave through the Maxwell and Poisson 
equations. This allows us to study the nonlinear interaction between arbitrarily large amplitude 
electromagnetic waves and a quantum plasma. We have used our system of nonlinear equations to 
study theoretically the parametric instabilities involving stimulated Raman scattering and modula- 
tional instabilities. A model for quasi-steady state propagating electromagnetic wavepackets is also 
derived, and which shows the possibility of localized solitary structures in the quantum plasma. Nu- 
merical simulations demonstrate the collapse and acceleration of the electrons in the nonlinear stage 
of the modulational instability, as well as the possibility of wake-field acceleration of the electrons 
to relativistic speeds by short laser pulses at nanometer length scales. The study has importance 
for the nonlinear interaction between a super-intense X-ray laser light and a solid-density plasma, 
where quantum effects are important. 

PACS numbers: 52.35.Mw,52.38.Hb,52.40.Db 

I. INTRODUCTION 

With the advent of the X-ray free-electron lasers [T] there are new possibilities to explore matter on atomic and 
single molecule levels. On these length scales, of the order of a few Angstrom, quantum effects play an important 
role in the dynamics of the electrons. Quantum effects have been measured experimentally both in the degenerate 
electron gas in metals and in warm dense matters [2J. It has also been found that quantum mechanical effects must 
be taken into account in intense laser-solid density plasma interaction experiments [3H5]. The interaction of large 
amplitude electromagnetic waves with the plasma can lead to various parametric instabilities [SHE] • At laser intensities 
around 10 19 W/cm and above, the nonlinearity associated with relativistic electron mass increase in short laser pulses 
plays a significant role. Furthermore, the relativistic ponderomotive force [9] of intense laser pulses produces density 
modifications. Thus, in a classical plasma, nonlinear effects associated with relativistic electron mass increase and 
relativistic ponderomotive force very important, since they provide the possibility of the modulational instability 
[TU1 HT] followed by a compression and localization of intense electromagnetic waves. In addition to the modulational 
instability, there are relativistic Raman forward and backward scattering instabilities [12H15] and the two-plasmon 
decay |16j instability that lead to strong collisionless heating of the plasma in the relativistic regime. The parametric 
instabilities of intense electromagnetic waves in magnetized plasmas have also been investigated [17H19) . 

However, for intense electromagnetic waves interacting with the plasma in the X-ray and 7-ray regimes, both 
relativistic and quantum effects must me taken into account on equal footing. Accordingly, in this paper, we present 
a simple nonlinear model, based on the Klein-Gordon (KG) equation coupled with the Maxwell equations that are 
capable of treating both the relativistic and quantum effects. Our work has applications in laboratories [2JI20], in 
quantum free electron laser systems [2TH23] , as well as in astrophysical settings [23] where white dwarf cores [25] and 
neutron stars |26j are strong sources of x-rays and 7-rays. 

The manuscript is organized as follows. In Sec. II, we present our mathematical model for the coupled KG and 
Maxwell equations, exhibiting nonlinear interactions between relativistic electrons and electromagnetic fields. Linear 
properties of the electrostatic and electromagnetic waves are discussed in Sec. III. Section IV shows hoe our governing 
equations lead to the wave equation that reveals the phenomena of relativistic self-focusing and relativistic self- induced 
transparency of electromagnetic waves. Section V is concerned with the theoretical and numerical investigations of 
the relativistic parametric instabilities in the quantum regime. Section VI deals with relativistic optical solitary 
waves. The nonlinear dynamics of interacting intense localized electromagnetic pulses, as well as the new phenomena 
of the formation of nonlinear Bernstein-Greene-Krushkal (BGK)-like modes and associated electron acceleration are 
described in Sec. VII. Section VIII contains a brief summary and conclusions. 
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II. MATHEMATICAL MODEL 



Historically, the Klein-Gordon equation (KGE) for an electron is obtained from the relativistic relation between the 
energy £ and the momentum p, viz. 

£ 2 = p 2 c 2 +m 2 e c\ (1) 

where c is the speed of light in vacuum and m e the electron mass. By the substitution £ —¥ ihd/dt and p — > —iKV 
in (1), where Ti is the Planck constant divided by 27r, we obtain the KGE for a free electron as 

h 2 ^-h 2 c 2 V 2 *P + myi; = Q, (2) 
where ip is the electron wave function. The free-particle KGE fulfills the continuity equation 

f + V.j = 0, (3) 



Pe = ~o o [ip -wr-iJ-^r ), (4) 



where 

2m P c 2 \ dt dt 
and 

i pit 

j e = — (r^- w). (5) 

We have multiplied the right-hand sides of Eqs. Q and ^ by the electron charge — e, so that p e can be interpreted 
as the electric charge density and j as the electric current density. Since p e is neither positive or negative definite, it 
cannot be interpreted as a probability density, however, it can be interpreted as a charge density which need not has 
a definite sign. 

We now wish to use the charge and current densities as sources for the self-consistent electromagnetic scalar and 
vector potentials </> and A for a quantum plasma. We, therefore, let ip represent an ensemble of the electrons. 
Introducing the electromagnetic potentials into the KGE, we make the usual substitutions ihd/dt —> ihd/dt + e(j> and 
—ihV — > —ihV + eA, obtaining 

W 2 ip - c 2 V 2 ^ - m 2 e c 4 ip = 0, (6) 
where we have defined the energy and momentum operators as 

W = ih^ t +e<j>, (7) 

and 

V = -ihV + eA, (8) 
respectively. The electric charge and current densities are now obtained as 



p e = -^ ± ^[rm+^(my}, (9) 

Zm e c : 



and 



respectively. We note that the charge and current densities obey the continuity equation 

t + v.j, = o. (11) 

The self-consistent vector and scalar potentials are obtained from the electromagnetic wave equations 



c 2 V x (V x A) + V-f =/ioOe, (12) 



dt 2 v ' dt 
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and 



r, 9A 1 / 



(13) 



where /io is the magnetic vacuum permeability and Eq is the electric permittivity in vacuum, and pi is the neutralizing 
positive charge density due to the ions. For immobile, singly charged ions, one can assume that pi — en , where n$ 
is the equilibrium ion number density. 



Using the Coulomb gauge V • A = 0, we obtain from Eqs. ( 12 1 and ( 13 ) 

8 2 A 



dt 2 



and 



c 2 V 2 A + V^- =Moc 2 j e , 



- — (pe + Pi 



(14) 



(15) 



respectively. Taking the divergences of both sides of Eq. ( 14 ), we have 



V 2 ^ = MoC 2 V-j e , 



so that Eq. ( 14 1 can be written as 



(PA 



c 2 V 2 A 



-yU C 2 V X (V X j e ). 



(16) 



(17) 



Equations (JsJ , (15) and (17) are our desired system that describes intense laser-plasma interactions in the quantum 
regime. 

The non-relativistic limit is obtained from Eq. |6j by substituting ip — exp(— im e c 2 t/h), and by using the 
condition \hd"i> /dt\ <C m e c 2 ^ , together with the normalization of ^ such that \E"f* — nQ is the electron number 
density at the equilibrium. In this limit, Eq. ([6]), yields the Schrodinger equation 



ih 



i 



dt 2m, 



■(-iKV + eA) 2 * + e0* = 0. 



(18) 



Here, and in what follows, we have used a simplified model and neglected the electron degeneracy pressure. The 
latter is important in dense matters where the electron degeneracy pressure appears due to the Pauli exclusion prin- 
ciple. For a non-relativistic plasma, the quantum statistical pressure has been introduced in a nonlinear Schrodinger 
model [27], but this has to be investigated for relativistic quantum plasmas. 



III. COLLECTIVE ELECTROSTATIC OSCILLATIONS AND FREE PARTICLES 

In the absence of the electromagnetic field (viz. A = 0), we still have electrostatic waves due to the charge 
separation between the electrons and ions. At short wavelengths, the quantum effects become important and give rise 
to dispersive effects in the electrostatic wave. At these wavelengths, there is an interplay between collective electron 
oscillations and free electron motion. When the wavelength is comparable to the Compton wavelength, the electrons 
become relativistic, and there are relativistic corrections to the dispersion relation for the electrostatic wave. 

In the derivation of the dispersion relation for relativistic electrons, it is convenient to first make the transformation 
ip — ip exp(— im e c 2 t/h), where the wave function ip obeys the wave equation 

q \ 2 — - 

ih— + m e c 2 + e0) i) + H 2 c 2 y 2 ^- mlc 4 ^ = 0, (19) 



dt 



and the electron charge density is 



^ ~d~t~^~dT i ~ I ) e ^l ■ (20) 
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We next linearize the system (19) by setting <fi = (j>i and ip = ipo + tpx, where <fii = <^>exp(iK • r — iilt) + complex 
conjugate, ipi = cxp(iK • r — iflt) + ip_ exp(— zK • r + iflt), and where |"0o| 2 — n o- Separating different Fourier 
modes, we obtain from (19) the dispersion relation for the electrostatic oscillations as £ = 1 + \e = 0, where £ is the 
dielectric constant and the electron susceptibility is 

= c 2 pe [4 m y~h 2 (n 2 -c 2 K 2 )} 
Xe h 2 (n 2 -c 2 K 2 ) 2 -4m 2 e c 4 n 2 ' 1 ' 

where uj pe = (n e 2 /s^m^) 1 ! 2 is the electron plasma frequency. We note that in the classical limit K — > 0, we have 
Xe = — w p e /^ 2 j while in the non-relativistic limit c — ¥ oo, we have Xe — ~u pe /{Q 2 — T^k 4 /Am 2 ). After some reordering 
of terms, the dispersion relation can be written as 

n 2 (n 2 - c 2 K 2 )(n 2 - c 2 k 2 - u 2 pe ) - Am 2 e c\n 2 - «%) = o. (22) 

In the classical limit h — > we obtain the Langmuir oscillations f2 = ui pei while in the limit c — > oo, we retain the 
non-relativistic result 



On the other hand, in the limit K 0, Eq. ([22j) yields two possibilities, one of which is the Langmuir oscillations 
at the plasma frequency, f2 = uj pe and the other one is oscillations with the frequency = 2m e c 2 /h. The latter 
corresponds to a negative energy state, which can be interpreted as positronic state. 
We note that there is a non-dimensional quantum parameter 

H = huj pe /m e c 2 (24) 



in Eq. (22 1 that determines the relative importance of the quantum effect. Typical values are H = 10~ 4 for the 
electron number density n e ~ 10 30 m~ 3 in solid density laser-plasma experiments and H — 0.007 may be representable 
of modern laser-compressed matter experiments with n e ~ 10 34 m~ 3 . This corresponds to oj pe = 8 x 10 16 s _1 and 
A e = 4 x 10" 9 m for H = 10~ 4 , and u pe = 5.4 x lO^s" 1 and A e = 5.5 x 10" n m for H = 0.007, where A e = c/uj pe 



is the electron skin depth. The non-relativistic result (23) is valid for electrostatic waves with wave numbers in the 
range 1 < K\ e < 1/H. For KX e < 1, the quantum corrections to uj « uj pe are different from (23) and turns the wave 
frequency slightly lower than uj pe . However, this effect is negligible for small values of H, and may be smaller than the 
degeneracy electron pressure effect, which is neglected here. On the other hand, the limit K X e > 1/H corresponds to 
relativistic particles with K > 1/Ac, where Ac = U/m e c « 3.9 x 10~ 13 m is the reduced Compton wavelength. For 
w„ e — > 0, we obtain the relativistic free particle dispersion relation 



n = T^- + \/^+c 2 K 2 , (25) 

where the upper sign (-) corresponds to the motion of a free electron and the lower sign (+) can be interpreted as the 
motion of a free positron. 



In Fig. 1, we have plotted the solutions of the dispersion relation (22) for H = 10 and H = 0.007. Both 
the electron plasma oscillations and the positronic states are shown. The electron plasma oscillations have a cutoff 
frequency lu — > cu pe when K — > 0, while the positronic states have a cutoff frequency lo — > 2m e c 2 /h, corresponding 
to uj/oj p e — > 1/H at K — > in Fig. 1. For the electron plasma oscillations, the increase in the wave frequency 
due to the quantum effect becomes noticeable approximately where KX e — 1/y/H, or K w (Airno/ 'as) 1 / 4 where 



m 



as = 47T£o?i /m e e 2 « 5.3 x 10 11 m is the Bohr radius. This corresponds to a wavelength of 2tt/K w 2.8 x 10 10 
for H = 10~ 4 and 2tt/K rs 5 x 10" 11 m for H = 0.007. 

The positronic states are associated with, for example, the Zitterbewegung effect [28, 29 , in which the interference 
between the positive and negative energy states are predicted to give oscillations on Compton wavelength scales in 
space. The Zitterbewegung effect is still debated and has not yet been observed in experiments. 



IV. NONLINEAR ELECTROMAGNETIC WAVE PROPAGATION AND SELF-INDUCED 

TRANSPARENCY 



It is well-known [30j that a large amplitude electromagnetic wave propagating in a classical plasma changes the 
dispersive properties of the plasma due to the relativistic mass increase of the electrons. We show here that the same 
effect occurs in our Klein-Gordon-Maxwell system. 
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FIG. 1: Dispersion curves (fi vs. K) for the electrostatic oscillations for H = 10 4 (solid curves), H = 0.007 (dashed curves), 
where H = hu) pe /m e c 2 . For K\ e > 1/H, the particle motion turns from weakly relativistic to ultra-relativistic. 



We consider for simplicity the propagation of a right-hand circularly polarized electromagnetic (CPEM) wave of 
the form A = Aq [x cos(koz — ujot) — y sin(fco2: — u>ot)] , where ujq is the wave frequency and fco the wavenumber. Due to 
the circular polarization, the oscillatory parts in the nonlinear term proportional to A 2 in the Klein-Gordon equation 
vanish. Assuming that tp depends only on time and not on space, and that = 0, we obtain from Eq. ([6]) 



dt 2 



+ m 2 e cS 2 A ip = 0, 



(26) 



where = yl + e 2 A 2 /m 2 c 2 can be interpreted as the relativistic gamma factor due to the electron mass increase 
in the CPEM wave field. Equation ( 26 ) has the solution 



tp = ip e-np(-im e c 2 ~f A t/h), 

where the constant tpo is determined by assuming the constant density p e 
Inserting (271 into Q with p e = — en Q , we obtain 



no 
J a' 



On the other hand, inserting (27) into (10) we have 

e 2 \M 2 



e 2 n 
lAm e 



which can be inserted into (14 1 to obtain 



-eno m 



Eq. j9j. 



(27) 



(28) 



(29) 



1 d 2 A 

d 2 ^ 



V 2 A = 



Equation ( 30 ) admits the nonlinear dispersion relation 



lA-m e 



A. 



(30) 



2 1.2 



pe 
1A : 



(31) 



which predicts a relativistic downshift of the CPEM wave frequency due to the relativistic electron mass increase in 
the CPEM wave field. Since the effective plasma frequency is decreased by a factor the model predicts the 

well-known self-induced transparency where the CPEM wave can propagate at frequencies below the electron plasma 
frequency. This is identical to the case of classical plasmas [3D] . 
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STIMULATED RAMAN SCATTERING AND MODULATION AL INSTABILITIES 



We now consider the instability of an intense CPEM wave in the quantum regime. In the presence of intense 
electromagnetic waves, we have the relativistic down-shift in the wave frequency given in (31 ), as well as the possibility 



of exciting electrostatic oscillations via the parametric instabilities. As an example, we will here consider stimulated 
Raman scattering instability, in which an intense electromagnetic wave decays into a daughter EM wave and an 
electron plasma wave. The two-plasmon decay instability, in which the CPEM wave decays into two electrostatic 
waves, will be treated elsewhere. 

It is convenient to first introduce the transformation tf> = •tpexp(—i"fA , m-eC 2 t/h), where ^a = yl + e 2 A^ j /m 2 c 2 and 
Aq is the amplitude of the EM carrier wave Aq. The wavefunction ip obeys the modified Klein-Gordon equation 



d 2 
^dt + ~ fArrieC + e <fi) ^ 



c 2 (~ 



ihV + eA) 2 -0 - m\c 



V = o, 



and the electron charge density is given by 



ihe 



2m P c 2 



dip dip 



dt 



dt 



1A 



(32) 



(33) 



Now, we linearize our system by introducing ip(r,t) — tpo + ipi (r, t ) (where ipo is assumed to be constant), A = 
Ao(r, t) + Ai(r, t), and <p(r,t) = <^>i(r, t). Using pi = eriQ into Eq. ( fl5| ), we note that the equilibrium quasi-neutrality 
requires that ipo is normalized such that IV'ol 2 — nail A- Using that Ao fulfills the plane wave equation (301, the 
linearized KGE (32 1, Poisson's equation (15) and the EM wave equation (17) then become 



d 2 ip 



- + c 2 V 2 ijji | + 2ih~i A m e c 2 ^$- + 2ihc 2 eA ■ Vipi 
) at 



+ ( 2j A m e c e<pi + ifte-gj- ) V'o - 2c e A • AiV>o = 0, 



(34) 



V 2 0! = 



ieh 



2eom e c 2 



dt 



dt 



and 



— ^Ai 

1A 



n 



i Vx{Vx[A (^Vi+lM'!)]} > 



(35) 



(36) 



respectively. We note that the term proportional to Ao ■ V^i in Eq. (34) gives rise to the two-plasmon decay, which 
we, however, do not consider here. 

We now introduce the Fourier representations ip = ip+ exp(— iSlt + zK • r) + rp_ exp(i£l£ — iK ■ r), <p\ — tj> exp(— ifM + 
zK • r)+ c.c, Ao = (1/2) Ao exp(— ujQt + ko • r)+ c.c, and Ai = [A + exp(— ioj + t + ik + ■ r) + A exp(— iuj-t + ik_ ■ r)] + 
c.c, where we introduced ui± = wo ± f2 and k± = ko ± K, and c.c. stands for complex conjugate. In one of the steps, 
we take the scalar product of both sides of the EM wave equation by Aq and use the fact that Aq ■ [k± x (k± x A )] = 
(k± x Ao) • (Aq x k±) = — |k± x Ao| 2 . Separating different Fourier modes and eliminating the Fourier coefficients, 
we find the nonlinear dispersion relation 



1 + 



(c 2 K 2 



Q 2 



-^IJia) 



Xe [4-/ 2 A m 2 c 4 -h z (n 2 ~c 2 K 2 )] 



c 2 e 2 |k+ x A | 2 



c 2 e 2 |k_ x A 



k 2 + D A (u + ,-k + ) k 2 _D A (oj-,k_) 



(37) 



where the electromagnetic sidebands are governed by Da((j->±, k±) = c 2 /? 2 . — to 2 . + uj 2 e /jA- The electric susceptibility 
in the presence of the laser field is given by 



h 2 (n 2 -c 2 K 2 



1A [K {V 2 - c 2 K 2 ) 2 - 4j%mlc 4 n 2 } 



(38) 
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After reordering of terms, the nonlinear dispersion relation (37) can be written as 

, ,2 



1 



pe 



rale 2 



D L (n,K) 



|k+ x A | 



k_ x A | 



k 2 + D A (u; + ,k + ) klD A (u>-,li-) 



= 0, 



where the electron plasma oscillations in the presence of the laser field are represented by 



J 2 

pe 
1A 



4-f A mlc 4 



n 2 



C 2 K 2 _ <^pe 
1A 



(39) 



(40) 



We note that Dt, = gives the dispersion relation for pure electrostatic oscillations in the presence of a large amplitude 
electromagnetic wave. 

In the classical limit fi — ¥ 0, we have \ e = —uj 2 e /fl 2 j A , and the nonlinear dispersion relation takes the form 



n 2 1A {c 2 k 2 



1 



Q 2 



47^777,2 c 2 



-pe ~lA""e 

which can be written in a more familiar form as 



e 2 |k + x A | 2 

k 2 D A (oJ+,k + ) + k 2 _D A (u_,k_) 



1 



c 2 K 2 



1 



pe 



A^ A m 2 c 2 



|k+ x A | 



|k_ x A | 



fciD A (w+,k+) k 2 _D A (u_,k_) 



= 0, 



(41) 



(42) 



with Dl = u) 2 e /~f A — fl 2 . These results can be compared with, for example, the dispersion relations obtained in Refs. 
[13), [T4"l [T6] for the relativistic case and in [6] for the non-relativistic case. 

To proceed with the numerical evaluation of the nonlinear dispersion relation, we choose a coordinate system 
such that the CPEM takes the form Aq = (x + iy)Ao and kg = fcoz, and, without loss of generality, we choose 



K = K {] z + K ± y. Then, we have K 2 = K 2 + K\, lA = yf 1 + e 2 \A \ 2 /m 2 e c 2 , |k± x A | 2 = [2(fc ± ^||) 2 + K\] \A \ 2 , 
K l{ ) 2 + K\. We also us 

yJc 2 kl+u) 2 J~f A . 



and k 2 ± — (feo ± ^||) 2 + K\. We also use that the carrier wave Aq obeys the nonlinear dispersion relation ojq 
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FIG. 2: Growth rate (f2j/w pe vs. isTn and Kx) f° r stimulated Raman scattering in the presence of a large amplitude CPEM 
wave, for the amplitudes ao = 1, ao = 5, and do = 20 (left to right panels) where ao = e\Ao\/m e c for H = 10 -4 (top panels) 
and H = 0.007 (bottom panels). 



We now assume that the wave frequency is complex valued, fi = Qr + iOj, where Qr is the real frequency and 
O/ the growth rate, and solve numerically the dispersion relation ( [39] ) for Q. In Fig. 2, we have plotted the growth 
rate for stimulated Raman scattering instability as a function of the wavenumbers Kn and K±, for a few values of 

a = e\A \/m e c and H — huj pe /m e c 2 . For all cases in Fig. 2, we used k c/uj pe = 20, which corresponds to a wavelength 
of 1.25 x 10 _9 m for H = 10~ 4 and to 1.7 x 10~ n m for H = 0.007, which is in the X-ray regime. We observe that 
for H = 10~ 4 , there is a broad spectrum of unstable waves, in particular for ao = 5 and ao = 10. For H — 0.007, 
we observe a reduction in the spectrum of unstable waves and the growth rate (relative to the electron plasma 
frequency) is slightly reduced. This is due to the fact that the wavelength of the unstable electrostatic oscillation 
approaches the critical wavelength, where quantum dispersive effects become important compared to the plasma 
frequency oscillations. For H = 10~ 4 , this wavelength is A crit = 27r/(47rno/as) 1 / 4 w 5 x 10~ 10 m, corresponding to 



a critical wavenumber K cr u = 1.25 x 10 10 m -1 , and for H = 0.007, we have \ cr it = 2n/(4irnoaB) ~ 2.8 x 10~ n m, 
corresponding to K cr n = 2.25 x 10 11 ni _1 . Hence, for H — 0.007 we have fco ss K cr n, which leads to the reduction of 
the growth rate due to the quantum dispersion effect. 

Furthermore, it should be mentioned that we do not find Raman-type instabilities involving the positronic states 
in Fig. 1. This is consistent from the point of view of the conservation of charges, since the production of positrons 
would violate the conservation of electric charges. 




FIG. 3: The growth rate (f2j/w pe vs. Ku and K±) for the modulational instability in the presence of a large amplitude CPEM 
wave, for the amplitudes ao = 1 (left panel), ao — 5 (middle panel) and ao = 20 (right panel) where ao — e\Ao\/m e c. We used 
H = 0.007 for all cases. 

In addition to stimulated Raman scattering instabilities, we have also the modulational instability that dominates for 
the pump frequencies uiq < 2u; pe / > /%4, and the corresponding wavenumbers fco < (w pe /c) y^S/jA ■ The modulational 
instability usually occurs for small modulation wavenumbers, and saturates nonlinearly by the formation of relatively 
small localized structures/solitary waves. In the past, such nonlinear structures have been studied for classical plasmas 
in ID [3T] and 3D [35] • We have investigated the modulational instability for the CPEM dipole field with fco = 0, and 
have plotted the results in Fig. 3 for different amplitudes ao = 1, ao = 5 and ao = 20. We find that the growth rate 
is relatively insensitive to the quantum parameter H . We have used H = 0.007 in Fig. 3, but H = gives almost 
identical results. This is understandable since the modulational instability takes place on relatively large scales and 
the quantum effect is thus negligible. However, we will investigate the quantum effect on the relatively small scale 
nonlinear structures below. 



VI. RELATIVISTIC OPTICAL SOLITARY WAVES 



Here we illustrate the existence of large amplitude localized CPEM wave excitations at the quantum scale in our 
system. We restrict our investigation to one-space dimension, which has also been studied for classical plasmas [31] . 

Far away from the local excitation, one can assume that the dynamics of the plasma is non-relativistic. To shorten 
the algebraic steps, it is convenient first to introduce a new wave-function ^>(z,t) and the potential $ via the trans- 
formations ip(r, t) — ^(z, t) exp(~im e c 2 t/h) and <fi = $ — m e c 2 /e, and which satisfy the KGE 



d 

l dt 



dz 2 



7 >^ 4 * = 0, 



(43) 



where 7a = (1 + e 2 v4 2 /mgC 2 ) 1 / 2 . In this gauge, the wave function \& is non-oscillatory in time, and the new potential 
takes the value $ = e/m e c 2 , far away from the solitary wave where the plasma is at rest. The electron charge density 
is expressed as 



m e c z 



i'h 
~2 



\[r* \Jf 

dt dt 



e$|#| 



(44) 



We now study quasi-steady state structures propagating with a constant speed v , so that (f) = </>(£) and A 2 — A 2 (£), 
where £ = z — v$t and A 2 = |A| 2 . The CPEM wave vector potential is of the form A = A(£)[xcos(fco;z — uJot) — 
ysin(fcoz — toot)]. It is convenient to introduce the eikonal representation = P(£) exp[i0(£)], where P and 9 are 
real- valued. Then, the KGE (43) takes the form 



h 2 (c 2 -v 2 ) 
— 2ihe&v 



d 2 P 
dP 



-p 



2i 



.dP d9 



iP- 



iP- 



} d6 
d£ 



d£ d£ ■ - de 

(e 2 $ 2 - m 2 c 4 7A )P = 0. 



ihevQ—P 



(45) 
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Setting the imaginary part of Eq. (45) to zero, we obtain 



where 



d6 



at, 



(46) 



(47) 



The solution of Eq. (46) is P U = D = constant. Using the boundary conditions P = no, 4> — (hence $ = m e c /e) 
and d/d^ = at |£| = oo, we have D = —nohm e c 2 VQ. Hence, we obtain 

d6 uo™e7o / e$ rig 



p2 



(48) 



where we have denoted 



7o = 



(49) 



The electron charge density ( 44 ) now takes the form 



Pe = 



hv 



d0 

di 



e$ P 



(50) 



which, by using Eq. ( 48 1 , can be written as 



v 2 e$ P 2 
.JL _| 

c 2 m e c 2 no 



(51) 



and hence Poisson's equation (15), with p$ = eno, becomes 

d 2 $ en 7o ( e * p 



d? 



- 1 



e \ m e c n 



(52) 



On the other hand, by setting the real part of Eq. (45) to zero, we have 

2" 



h 2 c 2 



d 2 p 



p 



d6\ 
d£j 



d6 

2he<S>voP-rz + (e 2 $ 2 - m 2 e c 4 jl)P = 0, 



(53) 



which, by using Eq. ( 48 1 , can be written as 

d 2 P 

h 2 — 



2 2 4 



i 2 c 4 c 2 P 4 



7o 2 



P = 0. 



(54) 



Finally, inserting the ansatz A = A(l;)[5L cos (koz ~~ w oO — ysin(fco.z — t+>ot)] for the vector potential, together with 
the current 



-|VfA 



-P 2 A, 



(55) 



into Eq. (17), we obtain the EM wave equation 

d 2 A uj 2 



d^ 2 



A + 76M1 



P 2 
»o 



A = 0, 



(56) 



where we used fc = w o w o/ c j an d where A = (c 2 /o;~ e )(a; 2 — a; 2 — c fc 2 )/(c — w 2 ) = w 2 /w 2 e — 7q is a nonlinear eigenvalue 
of the system that determines the wave frequency ujq. 
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FIG. 4: The spatial profiles of the vector potential, the electron charge density and the electrostatic potential (top to bottom 
panels), for H = 0.007 (left panels) and H = 0.002 (right panels) for standing solitary waves (vq = 0), and for A = —3.4. 
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FIG. 5: The spatial profiles of the vector potential, the electron charge density and the electrostatic potential (top to bottom 
panels), for H = 0.002, v = 0.0059c, and A = -3.4. 



The coupled system ( 52 1 , ( 54 1 and ( 56 1 describes the profile of electromagnetic solitary waves in a quantum plasma. 



It has the conserved quantity H = 0, where 



d e$ 
d£ m e c 2 



ll 



e 2 $ 2 P 2 
mlc 4 n 



27, 



, e<I> 

1 2 



n 1 



\ P 2 



dP 
1 



c 

l2~. 



p2 

n 



d eA 
dti, m e c 



li- 



lt 



p2 

n 



(57) 



The conserved quantity H can be used to check that the numerical scheme used to solve the nonlinear system (52 1, 
(54) and ( [56] ) produces correct results. 

In Fig. j^Twe have compared the present model with our previous results in Ref. [33j where we used a simplified 
model to describe the nonlinear interaction interaction of intense CPEM wave with a quantum plasma. We used the 
same parameters as in Fig. 2, of Ref. [33] to produce the profiles of the CPEM wave potential, the electron charge 
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density and the electrostatic potential. We observe that the present results are almost identical to our previous work 
[33] ■ For our sets of parameters, the quantum effect on the profiles of the solitary waves are small, and there is only 
a slight difference in the profiles of the electron density for the two values H — 0.002 and H — 0.007. For standing 
solitary wav es, su ch a s the one s in Fig. [4j the solutions are localized with exponentially decaying tails. By linearizing 
the system ( 52 ), (54) and (56) one can show that far away from the soliton, A decays as exp(V~ A£), while P and <fi 
are proportional to exp(A"£), where K is found from the dispersion relation 



hc 2 K 2 (c 2 K 2 



)-4 7o W e c> ^-u4) = 0. 



(58) 



For vq = (and 70 = 1), we see immediately that there exist only complex- valued K, which means that the quasi- 
stationary wave solutions decay exponentially far away from the solitary wave. However, in the classical limit h — > 0, 
we instead have the plasma wake oscillations given by real valued K — u pe /vQ. Hence, an electromagnetic pulse will 
create an oscillatory wake that extends far away from the EM pulse. In one-space dimension, there also exist special 
classes of propagating localized EM envelope solutions [M] [33]. In addition to the wake oscillations, we also have 
quantum oscillations in quantum plasmas. In Fig. 5, we show an example of a slowly moving envelope soliton, where 
small-scale oscillations in the charge density are clearly visible. 

We note that the cold fluid results can be retained in the classical limit fi — > 0. Then, Eq. (54) can be written as 



mfc? ~~ c 2 P 4 ' 7 2 



2 2 2 



By setting p e 



-en e , where n e is the electron number density, in Eq. (51) and solving for P, we obtain 

p2 



n 



1 



1 n 



75 \ n o 



1 



which can be inserted into ( 59 1 to obtain 



1 n e 



7o n o c 2 



n 



m e c 



1/2 



(59) 



(60) 



(61) 



which relates n e to $ and ja at a given speed vq. The relation (61) can also be obtained from the cold electron 



fluid model |34j and hence confirms the classical limit of the quantum model used here. If, furthermore, we assume 
standing waves such that Vq = 0, then we have from (61 1 



m e c 2 



1A- 



Solving for $ and inserting the result into Poisson's equation ([52]), we have 

1- 



A 



2^ 2 7A 



1A- 



"0 



(62) 



(63) 



where A e = c/ui pe is the electron skin depth. Finally, solving for P 2 and inserting the result into Eq. (56), we obtain 

,d 2 A loI . ( 2 d 2 -iA S\ A 



A 



de 



\A = 



de 



1 



1A 



(64) 



where we have used vq = and, therefore, fco = 0. Equation (64) is identical to the model of Marburger and Tooper 



[31j for the nonlinear optical standing wave in a classical cold fluid electron plasma. The relativistic mass increase is 



reflected by the ratio A/^a in the right-hand side of Eq. (64). The nonlinear electron density fluctuations, which are 
reflected in the term proportional to g^a/^ 2 in the right-hand side of Eq. (64) can often be neglected in the weakly 
relativistic case [55] . 

We note that our previous model [33j can be recovered in the weakly relativistic limit, in the following manner. 

Assuming that, to first order, we have a balance between the ponderomotive and electrostatic pressure so that 

.2 




1 + e(f>/m e c 2 « 7 a and that ja ~ 1, and v 2 <C c 2 . Accordingly, we have 



*2 2 
n c 



d 2 P 

He 



2m 2 c 4 1 



1A P = 0, 



(65) 
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1 , (66) 



Poisson's equation 

d 2 (j) _ eno / IaP 
d£ 2 s V «o 
and the CPEM wave equation 

We now make a simple change of variables "/aP = P- Then, we have, by neglecting terms containing djA/dS,, the 
model [53] 

h 2 c 2 + 2 ™'c 4 ( 1 + ~ 1a) P = 0, (68) 

and 



1 . (69) 



d£ 2 c 2 \7a' 

where the relativistic mass increase appears explicitly in the CPEM wave equation 



, xA=-^[- 1)A, (70) 

c 2 V7A"-o 



VII. NONLINEAR DYNAMICS OF INTERACTING ELECTROMAGNETIC WAVES IN QUANTUM 

PLASMAS 
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FIG. 6: The nonlinear stage of a modulational instability, showing localized electromagnetic wave-packets associated with 
depletions in the electron charge density, and positive electrostatic potentials. The electron pseudo-distribution function shows 
that the electrons have been accelerated and form BGK-like modes that travel away from the collapsed wavepackets. Parameters 
are H = 0.007, and initially a dipole field ao = 1 with ko = 0. 



In order to study the dynamics of the nonlinear interaction between intense CPEM waves and a quantum plasma, 
we have carried out numerical simulations of the Klein- Gordon-Maxwell system of equations. Wc have here restricted 
our study to one-space dimension, along the z direction in space, and written our governing equations in the form 

d \ 

ecj) V = W, (71) 



dt 

h^-+e<t>) W + h 2 c 2 
at 



d 2 *p 



j A m 2 e c 4 ^ = 0, 



(72) 
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1 d 2 A d 2 A 



c 2 dt 2 dz 2 



l^| 2 A, 



(73) 



and 



dz 2 



(74) 



We used a periodic simulation box in space, of length L x = 63 A e and used of the order 10 4 grid points to resolve 
the solution in space. It is important to resolve the relatively long electron skin depth scale as well as the shorter 
length scale associated with accelerated electrons with the momentum p = Tik and the associated wavelength A = 
27r/fc = 2i:h/p. Since we need at least two grid-points per wavelength to represent the solution, the required grid size 
is Ax < 7rh/p, which can be written Ax/X e < %Hm e c/p. For example, to represent the wave function of relativistic 
electrons with the momentum p = m e c, we need a spatial grid with Ax/X e < irH w 0.022 to represent the wave 
function for H = 0.007. The solution was advanced in time with the standard 4th-order Runge-Kutta scheme, using 
a timestep of order At = 10~ 4 Wp e . 
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FIG. 7: Attosecond laser pulse propagation into an underdense quantum plasma, at time t = 28.125 u)pe . Top to bottom panels 
show a) the electromagnetic vector potential of the laser pulse (the arrows show the propagation directions of the pulses), b) 
the electron charge density, c) the electrostatic potential, and d) the distribution of electrons in phase space in a 10-logarithmic 
color scale. Parameters are H — 0.007, amplitude ao = 1 and wavenumber ko = 20 A^ 1 . The laser pulses excite large amplitude 
oscillatory potential wakes behind them, as they penetrate the plasma slab. 
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FIG. 8: The same as in Fig. [7J at time t = 37.5uj pf }. Groups of electrons are accelerated to ultra-relativistic speeds by the large 
amplitude electrostatic wake neld. 
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We first study the growth and nonlinear saturation of the modulational instability, which is relevant for dense 
matters where the plasma is overdense or close to overdense. As initial conditions, we used a CPEM pump wave 
of the form A = A [x cos(fc z) — y sm (^o z )] with k — and A /m e c — 1. A small amplitude noise (random 
numbers) was added to A in order to seed any instability in the system. As initial conditions for the wave function, 
we used ip ~ \/noJ^A and W — m e c 2 ^/nojA, corresponding to a pure electronic state at equilibrium. The initially 



homogeneous electron density was set to n e = n = 10 m, corresponding to H = 0.007 [cf. Eq. ( 24 )] . In this 
situation, the electromagnetic wave is unstable due to the modulational instability, which has instability for small 
wavenumbers, as shown in Fig. 3. In the nonlinear stage, the EM waves self-focus into localized wavepackets similar 
to the ones in Fig. 4. Figure [6] depicts the late stage of the modulational instability. The collapse of the CPEM wave 
packet leads to relativistically strong ponderomotive potentials that accelerate the electrons to relativistic speeds. The 
relativistic electrons are associated with small-scale spatial oscillations in the wave function, where the wavelength 
is comparable to or even smaller than the Compton wavelength. We see in Fig. [6] that the CPEM wave envelope 
has been focused into localized wavepackets, associated with depletions in the electron density and positive localized 
electrostatic potentials. In order to study the distribution of electrons both in space and momentum space, we have 
performed a Fourier transform of the wavefunction ip using a moving window technique (using a Hann window) in 
space. The width of the window has been tuned so that it provides a good resolution both in space and in momentum 
space. The resulting spatial spectrogram gives a representation of the distribution of the electrons both in space and 
in momentum space; see Fig. |6][d), where the color indicates the density of electrons in phase space. In Fig. 6(d) the 
horizontal axis shows the spatial dependence and the vertical axis shows the momentum dependence via the relation 
p = hk between the momentum p and the wavenumber k. In this figure, it is clear that in the collapse stage of the 
solitary waves, bunches of electrons are accelerated to relativistic speeds and form self-trapped, Bcrnstcin-Grecne- 
Kruskal (BGK)-like modes that propagate away from the collapsed electromagnetic wavepackets. 

Next, we investigate a scenario of the short EM pulse propagation and the wake-field generation in a quantum 
plasma. This concept is traditionally used for the electron acceleration in classical plasmas [37l [38] . The numerical 
results are displayed in Figs. [7] and [8] Here, two atto-second pulses are injected from each side of a plasma slab 
and are allowed to collide at the center of the slab. As initial conditions, we used a CPEM pump wave of the form 
A = A (2;)[xcos(fcoz) — ysm(k z)] with fc = 20 A^ 1 and envelopes of the form A {z)/m e c — exp(— (z/X e ± 30) 2 ) 
propagating into the plasma slab. The plasma slab is initially centered between z — ±15A e with equal electrons with 
the number densities no, where the electron wave function was set to ip = y/n^ and W = m e c 2 y/n^. After a time 
t = 28.125 uipj-, we see in Fig. [7|b) and (c) that the large amplitude CPEM pulses excite plasma wake oscillations 
associated with large-amplitude positive potentials, and with an approximate wavelength of 2irc/u) pe , corresponding 
to a leading wavenumber of co pe /c. The positive potentials of the plasma wake oscillations are starting to capture 
populations of the electrons at edges of the plasma slab, at x ~ ±15A e . A high-frequency diffraction pattern is formed 
in the electron density, as faster electrons overtake slower electrons. Later, at t = 37.5 w^ 1 , the two laser pulses 
have collided and passed through each other. The trapped electrons have been further accelerated to ultra-relativistic 
speeds, as seen in Fig. pld) at x ~ ±5A e ., where the fastest electrons have reached a momentum of ~ 5m e c. 



VIII. SUMMARY AND CONCLUSIONS 



In this paper, we have developed a relativistic model for the interaction between intense electromagnetic waves 
and a quantum plasma. Our nonlinear model is based on the coupled Klein-Gordon and Maxwell equations for the 
relativistic electron momentum and the electromagnetic fields. In our fully relativistic model, the electron current 
and charge densities arc calculated sclf-consistcntly from the KGE, and they enter as sources for the nonlinear EM 
and electrostatic waves in the Maxwell equation. The KG-Maxwell system of equations has been used to derive 
the linear dispersion relation for the electrostatic and electromagnetic waves, as well as for investigating stimulated 
Raman scattering and modulational instabilities in the presence of relativistically intense CPEM waves. In the linear 
regime, the general dispersion relation for the electrostatic waves exhibits the quantum effect associated with the 
overlapping wave function. At long-wave-lengths, we have the dispersive Langmuir waves with frequencies close 
to the electron plasma frequency, while at shorter- wavelengths, we have the oscillation frequency of free electrons. 
At wavelengths comparable to or larger than the Compton wavelength, the electron motion is fully relativistic. 
In the nonlinear regime, we have demonstrated the existence of fully relativistic stimulated Raman scattering and 
modulational instabilities. While the Raman amplification is of much interest for generating a coherent electromagnetic 
radiation, the modulational instability gives rise to the localization and collapse of the CPEM waves into localized 
solitary EM wavepackets. Indeed, numerical simulations of the coupled KG-Maxwell equations reveal the collapse 
and acceleration of the electrons in the nonlinear stage of the modulational instability, as well as the possibility of 
wake-field acceleration of the electrons to relativistic speeds by short laser pulses at nanometer scales. In conclusion, 
we stress that the present investigation of nonlinear effects dealing with intense EM wave interactions with quantum 
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plasmas is relevant for the compression of X-ray free-electron laser pulses to attosecond duration 00] , as well as 
to the understanding of localized intense X-ray and 7-ray bursts that emanate from compact astrophysical objects 

nana- 
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